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The pseudogap state of high temperature superconductors is a profound mystery. It has tantalizing evidence 
of a number of broken symmetry states, not necessarily conventional charge and spin density waves. Here we 
explore a class of more exotic density wave states characterized by topological properties observed in recently 
discovered topological insulators. We suggest that these rich topological density wave states deserve closer 
attention in not only high temperature superconductors but in other correlated electron states. 



I. INTRODUCTION 

In a paper in 2000 Nayali^ provided an elegant classifica- 
tion of density wave states of non-zero angular momentum. 
The surprise is that given the roster of multitude of such states, 
so few are experimentally observed. Of these, the angular mo- 
mentum £ — 2, spin-singlet has taken on a special significance 
in the context of pseudogaps in cuprate high temperature su- 
perconductors,El It breaks translational symmetry, giving rise 
to a momentum dependent (1^2 _y2 (DDW) gap, without mod- 
ulating charge or spin, but alternating circulating charge cur- 
rents from a plaquette to plaquette much like an antiferromag- 
net. In its pristine form, in the half filled limit, that is, for 
one electron per site, the Fermi surface of DDW consists of 
four Dirac points and is therefore a semimetal. This broken 
symmetry state has inspired much effort in characterizing the 
pseudogap as a phase with an order parameter distinct from 
fluctuating superconducting order parameter. 

Presently, it appears from many experiments that the pseu- 
dogap may be susceptible to a host of possible competing or- 
ders. Thus it is important and interesting to explore an order 
parameter closely related to the singlet DDW, which retains 
many of its primary signatures such as the broken translational 
symmetry or a particle-hole condensate of higher angular mo- 
mentum. In particular we consider a density wave of non- 
zero angular momentum of mixed singlet and triplet variety 
such that in the half-filled limit, it is a gapped insulator Un- 
like the semimetallic DDW, it has a non-vanishing quantized 
spin Hall effect for a range of values of the chemical poten- 
tial. This is in fact a topological Mott insulator-' because it 
is the electron-electron interaction that is necessary for it to 
be realized. Further addition of charge caiTiers, doping, leads 
to Lifshitz transitions destroying the quantization but not the 
very existence of the spin Hall effect. 

It is remarkable that such an unconventional broken sym- 
metry, possibly relevant to high temperature superconductors, 
belongs to the same class of currently discussed novel state 
of matter known as topological insulators; in fact, our work is 
to some extent motivated by these recent developments.!^ We 
wish to emphasize that the undoped parent compounds of high 
temperature superconductors are proven to be antifeiTomag- 
nets with sizeable moments and the spin density wave trans- 
forms according to ^ = 0.^ The proposed topological density 
wave should therefore be relevant at larger doping that per- 



haps originates from a nearby insulating state. In no way is 
this different from the original suggestion of DDW. 

It has been known that triplet iadx2_y2 order parameter 
corresponds to staggered circulating spin currents around a 
square plaquette.* wherein the oppositely aligned spins circu- 
late in opposite directions, as shown in Fig. [T] This reminds 
us of topological band insulators where oppositely aligned 
edge-spins travel in opposite directions. However, there is no 
topological protection because the bulk is not gapped, but is 
a semimetal instead. A more interesting case is the order pa- 
rameter {icrdrc2_y2 + dxy), where cr = ±1 for up and down 
spins, with the quantization axis along z. Such a state not 
only satisfies time reversal invariance but is also fully gapped, 
analogous to time reversal invariant band insulators discov- 
ered recently. Singlet chiral {idx2_y2 +dxy) density wave that 
breaks macroscopic time reversal symmetry was employed to 
deduce possible polar Kerr effect and anomalous Nemst ef- 
fecP in the pseudogap phase of the cuprates. Another topo- 
logical state with a different symmetry of the order parameter 
was discussed in Ref. 8 

As to topological properties of superfluids, we refer the 
reader to the book by Volovik.^ Superconductors are particle- 
paiticle condensates, and, as such, the orbital wave function 
constrains the spin wave function because of the exchange 
symmetry. What we are discussing here are particle-hole con- 
densates, and there is no exchange requirement between a par- 
ticle and a hole. Thus, orbital wave function cannot constrain 
the spin wave function. Thus an orbital singlet can come in 
both spin singlet and triplet varieties. 

The plan of the paper is as follows: Section II is divided 
into three parts. Part A discusses the topological aspects in 
the absence of magnetic field, while Part B contains results 
for a perpendicular magnetic field. The Part C consists of a 
thorough discussion of the bulk-edge correspondence that fol- 
lows from topological considerations. In section III we dis- 
cuss Fermi surface reconstruction via a Lifshitz transition as 
the system is doped. In section IV possible experimental de- 
tection schemes are suggested. The symmetry of the order 
parameter that we have introduced is such that the necessary 
experimental techniques are more subtle than the detection of 
more common broken symmetries, such as spin or charge den- 
sity waves. 
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FIG. 1. (Color online) Triplet iad^2_y2 density wave in the absence 
of an external magnetic field. The current pattern of each spin species 
on an elementary plaquette is shown. The state is a semimetal. On 
the other hand iad^2_y2 + d^y can be fully gapped for a range of 
chemical potential. An example is shown in Fig. [2] 



with a generic set of band parameters, 

Efc = eife + £2/0 (6) 
eifc ~ —2t{coskx + cosfcj^), e2k = 4i' cos fc^. cos fc^ . (7) 

We may choose t = 0.15 eV, renormahzed by about a fac- 
tor of 2 from band calculations and t' ~ 0.3t, and Wq ^ 
—Aq ^ t ^ J, where J is the antiferromagnetic exchange 
constant in high temperature superconductors, for the purpose 
of illustration. Each of the two 2x2 blocks can be written 
in terms of two component spinors, ijjk,cr = {ck,a, Ck+Q^a)'^, 
a = ±1 = (t, i); for example, for the up spin block we have 

k 

(8) 

The eigenvalues (± refers to the upper and the lower bands 
respectively) 



II. ORDER PARAMETER TOPOLOGY 
A. Zero external magnetic field 

The order parameter that we consider is 

(4+Q,.Cfe.-'> = (<i>^(fc)r^)..,, (1) 

where cj, „{ck.a) is the Fermion creation(annihilation) opera- 
tor with momentum k and spin component ct; = 0, • • • 3, ti, 
T2, and Ta are the standard Pauli matrices and tq = 1. The 
nesting vector Q = [n / a^ir / a). We choose the components 
of the order parameter to be 

$^(/c) oc i^^{coskx — cos ky) = iWk (2) 

^°{k) (X Aq sin kx sin ky = Ak- (3) 

and the remaining components are set to zero. The right hand 
side is written in terms of the gap parameters and the con- 
version involves suitable coupling constants, which we do not 
need to specify in a non-selfconsistent Hartree-Fock theory. 
The lattice spacing a is set to unity. 

In the absence of an external magnetic field, the triplet d±id 
Hamiltonian is 



nd±^d- ^^N = J2'^lAk'i'k, 



(4) 



where the summation is over the reduced Brilloin Zone (RBZ) 

by ky ± kx = 



bounded by ky zL k^ — ±7r, and the spinor, ^j^., is defined 



as (4,t' 4+Q,r' 4,;' The chemical potential is sub- 

tracted for convenience, N being the number of particles.The 
matrix Ak is 



A. 



( efc-/i Ak + iWu 

Afe - iWk ek+Q - fj. 

ek-n Ak-tWk 

V Ak+iWk ek+Q-li 



(5) 



Afe,± = e2k — fJ- i Ek, Ek 



Ik 



(9) 



are plotted in Fig. [2] Since up and down spin components 
are decoupled, the Chern number for each component can be 
computed separately. After diagonalizing the Hamiltonian, 
we can obtain the eigenvectors 



-icrek/2^T 



where 



2^ eJ' 
1 eifc N 

Wk 

arctan(— — ) + 7r0(— Aj.). 
Ak 



(10) 

(11) 
(12) 

(13) 



To compute the Berry phase of the eigenstates, we define 
the Berry curvature, as 



(14) 



Substituting the eigenstates into the above equation, the 
Berry curvature can be written as 



^a,± = «Vfc X [(4 - 4)Vfc(«f^Y)]- 



(15) 



Since u±, w±, and 6k only depend on k^. and ky, only the z 
component, ricr,±, is non-zero, which is given by 



ilk ^ 

Ek 


d(h_ 
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dky 


dky 


Ak 


Wk 




dk^ 


dk^ 


dk^ 


dky 


dWk 
dky 


deik 
dky 



c—y- 



(16) 



From the above determinant, we can see that the Berry cur- 
vature will be zero if one of A^ and Wk is zero, so we need 



a mixing of (1^2 _y2 and d^y to have a non-trivial topological 
invariant. 

If we define the unit vector = ha-/\ha\, where — 
(Afe, —aWk, El), the Berry curvature can be written as 



I ^ , dfia drier , 
More explicitly, the Chern numbers are 



(17) 



a,± 



RBZ 



2lT 



[ 'L]L (sin^ ky + sin^ k^ cos^ ky) 

J RBZ 27r bj^. 



±cr. 



(18) 



We can focus on the lower band as long as there is a gap be- 
tween the upper and the lower bands. Then, 



So, the quantized spin Hall conductance will be 



spin 



' h 2e ^P'" 



e 

2^ 



(22) 



The eigenstates, Y^a,± {k)), are also the eigenstates of and 
Sz with eigenvalues 5^ = | and 5*2 — Since the spin 
SU{2) is broken by the triplet DDW, one might wonder if 
the Goldstone modes not contained in the Hartree-Fock pic- 
ture may not ruin the quantization. If SU{2) is broken down 
to C/(l), then there is still a quantum number corresponding 
to, say Sz, which is transported by the edge currents in the 
system. More succinctly, as long as time-reversal symme- 
try is preserved, we will still have Kramers degeneracy in our 
Hartree-Fock state, and therefore the edge modes will remain 
protected. 



N = A^i^,_ + A^;,_ = 



(-l)-l 



(19) 
(20) 



irrespective of the dimensionful parameters. Note, however, 
that the Chern numbers vanish unless both Aq and Wq are 
non-vanishing. The quantization holds for a range of chemical 
potential /i, as can be seen from Fig [2] 



B. Non-zero magnetic field 

In an infinitesimal external magnetic field, H, there will 
be a spin flop transition in the absence of explicit spin-orbit 
coupling, as shown in Fig. [s] We can assume H = Hz and 
the spins quantized along the x direction without any loss of 
generality. Then the Hamiltonian now becomes 




FIG. 2. (Color online) Energy spectra, Xk,± + jJ., corresponding to 
{iad^2_y2 + dxy) density wave. Here, for illustration, we have cho- 
sen Wo = t and Ao = —t and the band parameters, as described 
in the text. The chemical potential, /i, anywhere within the spectral 
gap, the lower band is exactly a half-filled and the system is a Mott 
insulator, unlike the semimetallic DDW at half-filling. 



For the fully gapped case, there will be a quantized spin 
Hall conductance associated with the eigenstates. The ratio of 
the dimensions of the quantized spin Hall conductance to the 
quantized Hall conductance should be the same as the ratio of 
the spin to the charge carried by a particle, since in two dimen- 
sions for both quantities the scale dependence L''^^ cancels, 
that is. 



(23) 



As before, the summation is over the RBZ, and the spinor is 
the same. The matrix Ak is now 






~Ak + iWk 




where ^ — ek + (t ^^^^ = <^1- Although the spin up 
and down components are coupled, particles with momentum 
k and spin up only couple to holes with momentum k + Q 
and spin down, and vice versa. Therefore, by redefining the 





-Ak - iWk 




Afe + iWk 





spinor, *J = (4,t'4+Q4'4,;'4+Q,t)' 
can still be expressed as a block diagonal matrix: 'Hd±id = 

\E'JAj.\['j,. The Chern numbers for each subblocks, i = 
1, 2, can be calculated as before. Therefore, defining 77^ = H- 1 
or -1 for i = 1 or 2, we obtain j = [(ei + 77^7)^ 

A2jl/2_ 



the Hamiltonian 



and the Berry curvature 
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FIG. 3. (Color online) Spins are flopped perpendicular to the applied 
magnetic field H. Contrast with Fig.[T| 



where hi — (77^ A^, — Wfe, ei + ijij). Performing a surface 
integration of the Berry curvature we get 



i,± 



cPk. 



RBZ 

At 



RBZ 27r-Efe^i 



sin^ k,, + sin^ fcr cos^ fc„ 



(cos kx sin ky + sin kx cos ky) 



(25) 



The ± refers to the upper and the lower band respectively. The 
integral does not depend on the external field, nor on magni- 
tude of the parameters t, W{), and Ao. The Chern numbers 
are 



0, 



Once again the spin Hall conductance is quantized, but the 
charge quantum Hall effect vanishes. The flopped spins carry 
the same current as before. The corresponding spin Hall con- 
ductance, as long as the gap survives, is 



spin 



27r' 



(27) 



The eigenstates, |<&i.±(fc)), are the eigenstates of with 



eigenvalues S 



\, but not eigenstates of because of the 



mixing of up and down spins. 



C. Bulk-edge correspondence 

For the {iadx^-y^ + dxy) order the bulk-edge correspon- 
dence can be studied by open boundary condition in the x- 
direction but periodic boundary condition in the y-direction, 
that is, by cutting open the torus. The edge modes if they exist 
will reside on the ends of the cylinders. The cut then leads to 
a Hamiltonian 



ky,i,j 



(28) 



where the spinor is ^fj^fc^ = (ci^fe^^f Q^fc^+Trf, Cj,fc„;Q^fe^+^;)^, 
and Aij{ky) is a AN x AN matrix parametrized by the wave 
vector ky, which is given by 



/ Tij{ky) Sij^^{ky) 

t 



Aij{ky) 







Sj^^^iky) T,,iky + n) 



V 



Tij(ky) Sij_^{ky) 

Sl^,^{ky) T,Aky + ^)J 



(26) where Tij{ky) and Sij,„{ky) are N x N matrices: 



Tij{ky) 



Sij^(j(^ky^ 



( 



-ji — 2t cos ky 
-t + 2t' cos ky 




-t + 



2t' cos ky 
- 2t cos ky —t + 2t' cos ky 
2t' cos ky — fi — 2t cos ky — t + 2t' cos ky 



V 



Wn 



I —2 cos ky —1 

2 cos ky 



1 




-1 




1 

—2 cos ky 



-t + 2t' cos ky 
\ 



t + 2t' cos ky 
2t cos ky I 



V 



(-1)^-1 

(-1)^ (-l)^2cosfcy/ 



.An 



sin ky 



1 

-1 
-1 



V 



-1)^ / 



The corresponding one dimensional system with N sites de- 
pends on the band structure and the order parameters defined 
above. 



The eigenvalue spectra are shown in Fig [4] The spectra, 
degenerate for up and down spins, are plotted in the range 
< ky < TT (ky < can be obtained by reflection). To find 
the edge states we choose the chemical potential in the gap. 
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In Fig. |4j we put /i = — 0.075eV for the purpose of illustra- 
tion. There are two edge states with positive group velocity, 
one with up spin and the other with down spin. Let them be 
ipy.^ and ipy^i, respectively. There are also two edge modes 
with negative group velocity denoted as V'<,t V'<,i for up 
spin and down spin, respectively. By explicitly computing the 
support of each of these wave functions, we have verified that 
electrons in states ijjy.i and tp<c.t localized near the left 
edge of the system whereas those in states 4'<,l ■0>.t 
localized near the right edge. The localization length of these 
states is essentially a lattice spacing; an example is shown in 
FigEl 
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FIG. 4. (a) Spectrum of the triplet (d± jd)-density wave on a cylinder. 
Parameters are t = 0.15eV, t' = 0.3t, fi = -0.075el/, Wo = t, and 
Ao — —t. The subscripts L and R to the spins correspond to left and 
right modes, (b) The probability density for positive group velocity 
for L and R spins for a lattice of A*' = 100 sites. 




FIG. 5. The bulk spectra for fixed values of ky with the same param- 
eters, as in Fig.|4] 



m. FERMI POCKETS AND LIFSHITZ TRANSITION 



It is interesting to track the evolution of successive Lifshitz 
transitions as we change the parameters. At first, when we 
lower the chemical potential, four hole pockets will open up 
in the full Brillouin zone, as shown in Fig. [6] and the corre- 
sponding spin Hall effect will lose its quantization but not the 
effect itself. But in mean field theory this cannot continue 
indefinitely with the nodal or the antinodal gaps fixed. So 
the parameters Wq and Aq will also decrease and will lead 
to a further opening of two electron pockets in the full Bril- 
louin zone, as shown in Fig. |6] Ultimately, when the doping 
is increased further, the large Fermi surface will emerge as a 
further Lifshitz transition. There is good evidence that such 
Lifshitz transitions indeed occur in high temperature super- 
conductors. 




-3-2-10 1 2 3 



-3-2-10 1 2 3 



It is interesting to see how this spectra compare with the one 
where periodic boundary conditions are applied in both x and 
y directions. After diagonalizing the Hamiltonian, we plot the 
spectra for a fixed value of ky for all values of the energies. 
The results are shown in Fig.|5] which are essentially identical 
to Fig.|4] except that the edge states are missing. 



FIG. 6. (Color online) (Left) Region plot, /i = —0.16 eV. Here for 
illustration, we have chosen Wo = t and Aq = — t as before. The 
hole pockets open up. (Right) Wo — O.OSt and Aq — —0.5t illus- 
trating the opening of the electron pockets at (vr, 0) and symmetry 
related points with enlarged hole pockets. 
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IV. EXPERIMENTAL DETECTION 

While there are many speculations about the nature of the 
pseudogap, they largely fall into two categories: 1) it is a 
crossover between a Mott insulator and a Fermi liquid, with- 
out any sharp, coherent excitations, and 2) it reflects a broken 
symmetry, with quasiparticles due to reconstructed Fermi sur- 
face that, despite strong correlations in the system, can behave 
in many ways as weakly interacting particles. The resolution 
of this dichotomy will ultimately be settled by experiments, 
which, to date, have shown some support for both. In the ab- 
sence of a definitive evidence one way or the other, we have 
adopted the second perspective (to some extent motivated by 
recent quantum oscillation experiments) to see what conse- 
quences there may be of having a broken symmetry phase 
with sufficiently hidden order, in particular one that has strik- 
ing similarities to topological insulators. 

A prime characteristic of a broken symmetry is that deep 
in the broken symmetry phase, an effective mean-field, or a 
Hartree-Fock Hamiltonian, suffices in discussing the proper- 
ties of matter, and the symmetries alone determine the ex- 
citation spectra and the collective modes. It is only in the 
proximity of quantum critical points that such a description 
breaks down but that is not the subject of discussion here. 
Moreover, those properties that are determined by symmetries 
alone should be robust and can be understood in the weak cou- 
pling limit, simplifying our task of exploring correlated elec- 
tron system. 

The mixed triplet-singlet order parameters considered here 
is even more hidden than the corresponding singlet DDW. Not 
only do they not modulate charge or spin, but so long as spin- 
orbit coupling is absent, they are also invisible to elastic neu- 
tron scattering because there is no associated staggered mag- 
netic field, as in a singlet DDW. 

Inelastic neutron scattering can detect its signature in terms 
of a spin gap at low energies in the longitudinal susceptibil- 
ity and signatures in the transverse susceptibility of quasi- 
Goldstone modes, and even onset of a finite frequency res- 
onance mode. Recall that at any finite temperatures SU{2) 
symmetry cannot be spontaneously broken in two dimen- 
sions; interlayer coupling is necessary to stabilize it. Thus 
the scale of symmetry breaking must be considerably smaller 
than t ^ J, and the signature must be sought at higher ener- 
gies. It could be a challenge to disentangle the signal from 
inelastic spin density wave excitations. On the other hand 
since the quasiparticle excitations are essentially identical to 
the singlet DDW, the quantum oscillation properties will be 
similar,"' except perhaps those in a tilted field," which is cur- 
rently being explored. The essence of this order parameters is 
modulation of spin current and kinetic energy. So, it will re- 
quire probes that can detect higher order correlation functions, 
such as the two-magnon Raman scattering. In the presence of 
modest spin-orbit coupling, it may be possible to find small 
shifts of nuclear quadrupolar frequency (NQR). The modula- 



tion of the kinetic energy arising from the d^y component, in 
particular staggered modulation of t', may lead to anomalies 
in the propagation of ultrasound^lat a temperature where such 
an order is formed, presumably at the pseudogap temperature 
T* . The detection of the unique features of the proposed or- 
der parameter, the spin Hall effect and edge currents would be 
even more challenging. 

The effects of non-magnetic impurities on the mixed triplet- 
singlet phase studied here are rather subtle. We expect such 
disorder to couple only weakly to spin currents. Generically, 
disorder will couple differently to the iad,j.2_y2 and d^y com- 
ponents since each breaks a different symmetry. However, by 
breaking both the point group and lattice translation symme- 
tries, disorder can enable mixing with (generally incommen- 
surate) density wave states in other angular momentum chan- 
nels. For example, at the level of Landau theory, we expect 
terms in the free energy proportional to product of quadratic 
powers of the component order parameters, which would be 
proportional to the impurity concentration, thus inducing spin 
or charge density waves. So long as spin rotational symmetry 
is preserved in the normal state, the phase transition into the 
iad,j.2_y2 state can remain sharp. 

From the standpoint of topological order at zero tempera- 
ture, the effects of weak disorder are somewhat simpler. Since 
the density wave phase considered here is a gapped phase with 
topological order that is protected by time-reversal symmetry, 
it remains robust against weak non-magnetic disorder. Thus, 
the phase can still be described in terms of its topology at zero 
temperature, a feature which it shares with topological band 
insulators. 

Lastly, we remark that in the presence of magnetic impu- 
rities, the phase is not sharply defined - either as a broken 
symmetry or in terms of it's underlying topology. 

In terms of microscopic models beyond the phenemenol- 
ogy discussed here, it is almost certain that correlated hopping 
processes will play a key role,'^ Finally, since dx2_y2 and dxy 
are two distinct irreducible representations on a square lattice, 
generically they will each have their own transition tempera- 
tures, as dictated by Landau theory. The development of the 
dxy order parameter would be at a higher temperature com- 
pared to the triplet component which breaks SU (2) and there- 
fore requires interlayer coupling. Thus it follows that when 
applied to cuprates there must be two transitions in the pseu- 
dogap regime. Since the topological phase studied here arises 
from spontaneous symmetry breaking, it can support charged 
skyrmion textures in analogy with.''* The properties of such 
textures and their transport signatures shall be the topic of a 
forthcoming publication. 
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